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Abstract. The transmission through a stack of identical slabs that are separated 
by gaps with random widths is usually treated by calculating the average of 
the logarithm of the transmission probability. We show how to calculate the 
average of the transmission probability itself with the aid of a recurrence relation 
and derive analytical upper and lower bounds. The upper bound, when used as 
an approximation for the transmission probability, is unreasonably good and we 
conjecture that it is asymptotically exact. 



PACS numbers: 42.25.Dd 

1. Introduction 

We revisit a classical problem: The transmission through a linear array of many 
identical slabs (glass plates, plastic transparencies, or the like) with random separation, 
as depicted in Fig. 1. The transmission probability that Stokes derived in 1862 [1] on 
the basis of ray-optical arguments (thereby improving on an earlier attempt by Fresnel 
in 1821; see Refs. [2] and [3] for the history of the subject) is not correct because there 
are crucial interference effects that require a proper wave-optical treatment. Just that 
was given by Berry and Klein in 1997 [3] who found that the average of the logarithm 
of the transmission probability through N slabs is equal to N times the logarithm of 
the single-slab transmission probability 

(logr w ) =Nlog Tl . (1) 

Here, t y is the probability of transmission through a single slab and t n denotes the 
transmission probability for N slabs. Its implicit dependence on the random phases 
that originate in the random spacing of the slabs is averaged over, indicated by the 
(■ • •) notation. As emphasized in Ref. [3], the disorder is crucial; without it, most 
wavelength components would be transmitted, and the stack should then appear 
rather transparent, but this is not the case as a simple experiment with a stack of 
transparencies demonstrates [3, 4]. 
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Figure 1. A stack of JV identical slabs, each with single-slap transmission 
probability Tj. The stack as a whole has transmission probability r N , which 
depends on the phases that result from the random spacing of the slabs. We are 
interested in (t n ), the transmission probability of the stack averaged over the 
JV — 1 phases. 



It is indeed common to average logarithms because they are known to be "self 
averaging" [5], and the exact result (1) is truly remarkable. But one should realize what 
it tells us about the average transmission probability (t n ) itself. As a consequence of 
the inequality 

(logr*) <1o S {t n ) (2) 

the Berry-Klein relation (1) amounts to a lower bound on the average transmission 
probability, 

(^>>rf. (3) 

As we shall see below, this bound is not particularly tight because there is a very large 
range of individual t n values. In particular, we note that the ray-optics result [3] 

= r 1 +JV( 1 l-r 1 ) (4) 

is consistent with (3). 

It is the objective of the present contribution to report good wave-optics estimates 
for (t n ^ and closely related quantities. In particular, we will improve on the lower 
bound of (3) and supplement it with an upper bound. We observe that the upper 
bound, when used as an approximation for (t n ), is unreasonably good and seems to 

give us the exact asymptotic values of quantities such as (r N+1 } / (t n ) or (t^) 1 ^. At 
present, this coincidence of the upper bound with exact asymptotic values is a poorly 
understood mystery. 

2. Single slab: The transfer matrix 

For a wave of wavelength 2w/k, the wave functions to the left and to the right of the 
nth slab are 

^ ri « ht ) (x) = u n e ik( - x ~ Xn ~ ^ + v n e~ ik ( x ~ x " ~ £ ) , (5) 
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Figure 2. Amplitudes on both sides of the nth slab. The unitary scattering 
matrix S of (6) relates the incoming amplitudes « n -i and v n to the outgoing 
amplitudes u n and v n —i, whereas the transfer matrix T of (8) connects the 
amplitudes on the left with the amplitudes on the right. 



where x n is the position of the left edge and t is the thickness of the slab; see Fig. 2. The 
incoming amplitudes are related to the outgoing amplitudes by the unitary scattering 
matrix S 1 , 

ib' 



S 



U n -l 



(6) 



Vn-l I \ f, 

where the entries of S are restricted by 

\a\ 2 = \af =r 1 , \b\ 2 = \bf = 1 - t x , a*b' + b*a' = 0, (7) 

which account for the single-slab transmission probability t x and the unitary nature 
of S. The particular values of the complex phases of a, b, a', and b' are of secondary 
interest, but we note that we have a = a' = e lk ^ and b = b' = for a completely 
transparent, non-scattering slab, for which 2pn Ch \x) = ^n ieht \x). 

The transfer matrix T is used to express the amplitudes on the right in terms of 
the amplitudes on the left, 



T 



Un-l 
V n -1 



(8) 



The one-to-one relation between S and T implies that the transfer matrix is of the 
form 



where 



T 
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(9) 



(10) 



1 cosh(26>) + 1 ' 

and a, /3, f3' are phase factors that have fixed values which, however, are largely 
irrelevant for what follows. 

The transfer matrix for the gap of length L n between the nth slab and the (n+l)th 
slab is the diagonal phase matrix 

/ AkL n 



D(kL n ) = 









-ikL n 



(11) 
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Phase matrices of the same structure sandwich the central ^-dependent matrix in (9), 
so that we have 

/ cosh# sinh0 \ 
T = e ia D((3)t(6)D(f3') with t{6) = (12) 

V smhft coshft I 

as a more useful way of writing T. 

The product of two transfer matrices is another transfer matrix, whereby the 
relevant observation is the composition law 

t(0i) D(<f>) t(9 2 ) = D( 7 ) i(0 1&2 ) D{i) (13) 

with 0!&2 determined by 

cosh(26»i &2 ) = cosh(26»i)cosh(26» 2 ) 

+ cos(20)sinh(26>i)sinh(26> 2 ) (14) 

and the phases 7 and 7' by 

e 2i7 sinh(2<9i fc2 ) = sinh(20i) cosh(26> 2 ) 

+ cos(20) cosh(26»i) sinh(26> 2 ) 

+ isin(2</>) sinh(2# 2 ), 
e 2 ' 7 ' sinh(20 1&2 ) = cosh(26»i)sinh(26» 2 ) 

+ cos(20) sinh(20i) cosh(26> 2 ) 

+ isin(20)sinh(26»i). (15) 

Whereas (15) is of no consequence for the following considerations, the rule (14) is of 
central importance. 

3. Many slabs: A recurrence relation 

We now turn to the situation of Fig. 1, where we have N identical slabs separated by 
gaps L\, i 2 , . . . , Ln-i that are not controlled on the scale set by the wavelength 27r/fc. 
Therefore, we regard the phase factors e lkLn as random with a uniform distribution 
on the unit circle in the complex plane. 
The over-all transfer matrix 

T tot = TB{kL Y )TD(kL 2 )T ■ ■ ■ TD(kL N ^)T 
= e iNa D(p) ^ J] t(6) Dfo^J t{9) D(P') 

= c ia ^D(l3 tot )t(0 tot )D(l3> tot ) (16) 
is characterized by 6 to t which is obtained by repeated application of the composition 
rule (14), whereby the phases <j>„ = (3 + /?' + kL n have random values. Each 
experimental realization of the TV-slab stack of Fig. 1 has different values for these 
random phases, and the transmission probability 

2 

Tn = cosh(2# tot ) + 1 (17) 

varies from one experiment to the next. We need to average over the N — 1 random 
phases to find (t n ). 

Let us consider a somewhat more general question: What is the average value 
(/(cosh(20 t ot))) 01 a function of cosh(2# to t), and thus of a function of t n 1 When 
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Figure 3. Regarding the meaning of / n (C") in (18) and (19). The random phases 
4>i, <t>2, • ■ ■ , <t>n—i have already been averaged over, but the averaging over the 
phases <j> n , . . . , 4>n—1 is yet to be performed. 



the averaging is carried out successively, first averaging over <j>\ , then over ^2 , and so 
forth, finally over 4>n-i, we have an intermediate value f n (C') after averaging over 
the first n—1 random phases; see Fig. 3. Here C denotes the value of cosh(2# n ...jv) 
for the stack of slabs n to N with its dependence on the remaining phases <f> n , . . . , 
4>n-i- We then have fi(C) = /(C) for the value prior to any averaging, and (14) 
tells us that we get f n+ i(C) from f n {C) by means of 



f n+1 (C')= J ^/„(CC" + SS> cos <p) 



(2tt) 



with C = cosh(26<) = 1 , S = sinh(26>) = —y/l - t x , 



and S' = \JC' 2 - 1 with C" > 1 , (18) 

and the integration covers any convenient interval of 27r. Eventually this takes us to 

(f{cosh(29 tot ))) = f N (C) (19) 

when the recursive averaging over the N—1 random phases is completed. 

For illustration, we take /i(C') = C as a first example. The recurrence relation 
(18) yields f n (C) = C^C, so that 

( cosh(20 to t)) = cosh(26i) A ' (20) 

or 



^/r N ) = \ + \{2/r 1 -l) 1 



(21) 



when stated in terms of transmission probabilities. A second illustrating example is 
/ 1 (C") = C*' 2 -i, for which 



fn(C) 



fi(C), 
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(22) 



(23) 



follows. 

Taken together, (20) and (23) tell us that the normalized variance of cosh(2# to t) 
grows exponentially with the number of slabs, 



cosh(2fl t 



1 



(cosh(20 tot ))' 



cosh(20) 
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(24) 
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The values of cover a correspondingly large range, and so we understand why the 
two sides of (2) differ by much. 

This brings us to the much more important log r case of 

/i(C)=log-A-. (25) 



C' + l ' 



Here, 



f n (C) = (n - l)A(C) + A(C') , /^(C) = Nh{C) (26) 
is a manifestation of the "self-averaging" of the logarithm (not any logarithm though, 
but this particular one), and we get 

i°g .„u,J ^ i ; ) = Sn (C) - log n . (27) 



cosh(20 to t) + 1, 
This is the Berry-Klein result (1), of course. 

Finally, we turn to calculating (t n ). The first few f n (C')s are 

h(C) 1 



/2(C) = 
/3(C) = 



C' + l ' 
2 

C' + C 



^(C' + 1)(2C 2 + C'-1) ' 



(28) 



giving 



/ 2 (C) = £ 



(r 3 ) = / 3 (C) = 



2- Tl ' 
2 



(29) 



(C + 1)\/2C^T ^4/r x - 3 ' 
and it is frustratingly difficult to go beyond n = 3. But it is possible to evaluate the 
recurrence relation (18) numerically and so determine (j N ^ — f N (C). In passing, we 
note that ( T 2) rliy = ( T 2) an d ( T 3) ra y > ( T 3) ^ or < T i < 1; ra y optics fails for N > 2. 

For r x = 0.85, the outcome of such a computation is shown in the lin-log plot of 
Fig. 4 as the dotted curve 'a'. The crosses near the curve were obtained by a Monte 
Carlo calculation in which 400, 000 experiments were simulated with up to 200 slabs. 
The straight dashed line 'b' is the lower bound of (3). The solid lines are the upper 
and lower bounds discussed in the next section. Other values of r x result in plots with 
the same general features. 
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Figure 4. Average transmission probability for a stack of N identical slabs. 
For Tj = 0.85 and N < 200, the dotted curve 'a' shows the values of log^r^^ 
computed by a numerical evaluation of the recurrence relation (18), commencing 
with the small- ra functions of (28). The crosses that follow curve 'a' are values 
obtained by a Monte Carlo calculation that simulated 400, 000 experimental 
realizations. The two solid lines are the upper and lower bounds of (32) and 
(36), respectively. The dashed line 'b' is the lower bound (3) derived by Berry 
and Klein [3]. 



4. Many slabs: Upper and lower bounds 

Since / 2 (C") < l/VC^C, we have 

1 f dtp 



< 



(2tt) 
1 f dip 



2tt y/C + S(S'/C) cos <p 



^CKj J 2tt ^/C + Scoslp 

(27T) 

= -4=T(t 1 ), (30) 

where the second inequality recognizes that the integral in the first line is a 
monotonically increasing function of S'/C, so that the replacement S'/C — > 1 
increases its value. The integral defining T(r 1 ) is of elliptic kind and its value is 
less than 1 if C > 1, that is: T(r 1 ) < 1 if t x < 1. We conclude by induction that 

f n (C) < -^T^)"" 2 (31) 

holds for n > 2. The upper bound 

(r N ) < (T 2 )r( Tl ) N - 2 (32) 
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Figure 5. Upper bound T(r 1 ) and lower bound A(r 1 ) on ((T Ar }/(r 2 )) 1 /( JV 2 ' as 
functions of t 1 . The dashed straight line is the lower bound on (t^) 1 ^ of (3). 



then follows. The ray-optics result (4) is inconsistent with this upper bound. Figure 5 
shows T(r 1 ) as a function of t 1 . 

We derive a lower bound by first observing that 



with 



. h{C) 
m - Men 



for n > 2 



y/n ~ rf/4 for r x < 2 



V2 
V2 



(33) 



(34) 



(35) 



(36) 



and then inferring by induction that 

/JCO^coa^)"- 2 

holds for n > 2. The lower bound 

then follows. The plot of A(r n ) as a function of t x in Fig. 5 shows that A(tt) > r x for 
< T]_ < 1 and, therefore, this lower bound is more stringent than (3), but it is not 
tight either. We are certain, however, that (t n ) is bounded exponentially both from 
above and from below. 

Figure 6 illustrates the two bounds 

l/(N-2) 



A(7i)< 



for 



0.85. 



The values for curve 



<T(rJ 



(37) 



are obtained by the numerical evaluation of 



the recurrence relation (18). Clearly all values are well within the two bounds, the 
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Figure 6. Values of ((r Ar >/(r 2 )) 1 /( Ar - 2 ) for t x = 0.85. The bounds of (37) are the 
two horizontal dashed lines. Curve 'a' shows the actual values. The extrapolation 
explained in the context of (38) and (39) gives curve 'b'. 



horizontal dashed lines. This figure, and analogous plots for other values of t x , suggest 
that 

l/(AT-2) 

— ►T(t 1 ) asTV^oo. (38) 

The corresponding observation in Fig. 4 is that, for sufficiently large N, line 'a' there 
is parallel to the solid line for the upper bound. At present, (38) is no more than a 
conjecture that is supported by a body of numerical evidence. 

Some of the evidence is curve 'b' in Fig. 6. Its values are obtained by an 
extrapolation that assumes that 

l/(JV-2) 

ksA- B/N (39) 

for large TV with A and B slowly varying with N. For two consecutive N values of 
curve 'a' we can get an estimate of A and B, and curve 'b' represents the successive 
values of A thus extrapolated. The rapid and consistent approach of 'b' to the 
horizontal line of the upper limit feeds the expectation that the conjecture (38) could 
be true. We leave the matter at that. 

5. Summary 

We established the recurrence relation (18) that facilitates the calculation of the 
average value (f(r N )) of any function of t n , the transmission probability through 
the stack of N identical slabs with random gaps between them. We observed that 
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the individual values of t n are spread over a large range and, therefore, (t n ) exceeds 

e (logr JV > = T N by much ^ 

Further, we derived strict upper and lower bounds on (t n ^, both bounds being 
exponential functions of N. The ray-optics prediction for (t n ) is consistent with the 
lower bound but not with the upper bound. The upper bound, when used as an 
approximation for (t n ), is of much better accuracy than its derivation suggests and, 
based on numerical evidence, we conjecture that it is asymptotically exact. 
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